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Abstract 

By dressing Polyakov lines with appropriate functionals of the gauge fields, 
we construct observables describing a fundamental representation static quark- 
antiquark pair in the singlet, adjoint and average channels of SU(N) pure gauge 
theory. Each of the potentials represents a gauge invariant eigenvalue of the 
Hamiltonian. Numerical simulations are performed for SU(2) in 2+1 dimensions. 
The adjoint channel is found to be repulsive at small and confining at large 
separations, suggesting the existence of a metastable (N 2 — l)-plet bound state. 
For small distances and temperatures above the deconfinement transition, the 
leading order perturbative prediction for the ratio of singlet and adjoint potentials 
is reproduced by the lattice data. 



The potential of two static color charges in the fundamental representation separated 
by a distance r is a quantity of fundamental importance for the study of the confining 
force in SU(N) Yang-Mills theory and QCD, as well as the phenomenology of heavy 
quark systems. A quark-antiquark pair at the same point can be in either a singlet or 
an adjoint state, acccording to the irreducible representations of the tensor product 

N ® N = I + (N 2 — 1). (1) 

In perturbation theory, one fixes a gauge and a separate description for the two channels 
is possible also for small separations between the charges. To leading order one finds 
the relation (0(g 2 ) contributions cancel) |T], |2j 

Vi(r) = -(N 2 - l)Vjv»-i(r) + 0(g 4 ). (2) 

Because of the growing interaction strength, perturbation theory can be trusted only at 
short distances, while non-perturbative methods such as lattice simulations are required 
in the regime of linear confinement. 

On the other hand, in lattice gauge theory one is interested in non-perturbatively 
gauge invariant quantities onlyQ. It has so far not been possible to construct such 
operators for the two channels separately. At zero temperature, the standard lattice 
operator to compute the static potential is the Wilson loop projecting onto the singlet 
channel. Good agreement between an integration of perturbatively truncated renomal- 
ization group equations and continuum extrapolated lattice results have been reported 
for separations up to « 0.3fm [Q. At finite temperature, the only known gauge in- 
variant operator related to the potential is the Polyakov loop, which in turn yields a 
weighted average of the singlet and adjoint channels |5|, 0]. Thus, at asymptotically low 
temperatures, there are separate gauge invariant operators for the singlet and average 
potential. Together with the demonstrated validity of perturbation theory at short 
distances, this suggests that the octet channel exists also non-perturbatively. 

Recently it has been shown that, at the expense of working with non-local functionals, 
it is possible to arrive at a gauge invariant and non-perturbative description of two- 
point functions for charge carriers, whose exponential decay is governed by eigenvalues 
of the Hamiltonian ||. This has been used to define and compute a parton mass for 
the gluon, which can be related to level splittings of static mesons 0. 

In the present letter, we apply this formalism to Polyakov loop correlators in order 
to establish the non-perturbative existence and gauge invariance of the static potential 

1 Attempts were made to address this relation by lattice simulations in a fixed Landau gauge || . 
However, besides problems with uniqueness, lattice Landau gauge violates positivity and it is not clear 
whether the results are associated with physical quantum mechanical states. 
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in the singlet, adjoint and average channels. After briefly summarizing the symme- 
tries of a quark-antiquark pair, non-local eigenf unctions of the spatial lattice Laplace 
operator are used to construct a dressed Wilson line, which is gauge invariant up to a 
global residual symmetry of the Laplacian. Two such dressed Wilson lines separated 
by r represent an alternative operator to extract the zero temperature singlet poten- 
tial, producing identical results to those obtained from the Wilson loop. In the finite 
temperature case, the dressed Polyakov loop allows for separately gauge invariant ex- 
pressions for the average and the singlet potentials, and the adjoint potential may be 
computed from the appropriate difference between those operators. The transfer ma- 
trix formalism is used to show the gauge invariance of the potentials, and a numerical 
test in 2+1 d SU(2) Yang-Mills theory demonstrates the feasibility of simulating the 
operators in question. 

We consider SU(N) Yang-Mills theory in d + 1 dimensions with Wilson action on a 
L d ■ N t lattice. For a given Euclidean time, we are interested in the state with a static 
quark at y and its charge conjugate at x, ^ a (x)^ j g(y). The singlet and adjoint states 
of the pair are obtained by means of the projection operators |], § 

Pt = ^l®l-|f a ®T a , 
N 2 - 1 2 

Pad = — ^-l®l + -^T a <g)T<\ (3) 

where f a = -T a * and Pfa = P 1M , P^ = P ad P x = 0, P ad + P x = 1 <g> 1 . Applying 
these operators to the quark-antiquark pair, one finds the tensor decomposition 

(PaV5(x)V>(y)) a/3 = ^«/A( x )^7(y) 

P ad ^(x)^(y)) a/3 = ^a(x)^/j(y) - ^5 Q/3 ^ 7 (x)?/; 7 (y) 

= -2f a > 7 (x)T 7 ^^(y). (4) 



( 



The projection operators act on the space of representation matrices, but are space-time 
independent. Therefore, under local gauge transformations, the expressions Eq. (f|) 
behave as singlet and (A^ 2 — l)-plet for x = y only. 

The propagation of a static quark in Euclidean time is described by the (untraced) 
temporal Wilson line M, 



t-i 



V>a(x,t + t) = L(x) a/3 ^ /3 (x,to), £(x) = JJ [/o(x,*')> ( 5 ) 

t'=t 

while that of a quark-antiquark pair is given by the matrix correlation function Cq q (r, t) 
of two Wilson lines separated by r = |x — y|. According to the tensor decomposition 
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this can be written as a sum of singlet and adjoint channels, 

Cg q (r,t) = (Lt(x)L(y)) = e- tv ^P 1 + e- tv ^P ad . (6) 
Applying the projection operators one solves for the respective channels to find 

e -tVL(r) = Tr(P 1 Q (? )/TrP 1 = l(TrLt(x)L(y)), 
e -tv ad (r) = T r(P ad Cy/Tr(P ad ) 

= jy^I^^WT^ly)) ~ N(N l _ r) (TrL\x)L(y)). (7) 

Again, these expressions are only gauge invariant for x = y. This lack of gauge 
invariance is a deficiency of the operators: so far they do not reflect the fact that 
between the qq pair a string of flux is formed compensating for the transformation 
of the quarks. However, the singlet expression can be interpreted as arising from the 
gauge invariant Wilson loop in a particular gauge |9], |I|]. For example, evaluating a 
Wilson loop in the (xj, t)-plane in an axial gauge, where all links in the i-direction are 
fixed to the identity, one obtains the chain of equations 

(W(r,t)) = (W(r,t)6[Ui(x),l]) 

= (L^(x)C^(x,y)2^(y)E^ a (x,y)«5[^(x), 1]) 

= (Li /3 (x)^ 7 L 75 (y)^ Q ) = (TrLt(x)L(y)), (8) 

where Z7j(x, y) denotes the straight line of links in the z-direction connecting x and 
y. It represents the physical flux tube formed between the charges, and the whole 
quark- ant iquark-glue system is in a singlet state. At T = 0, this is the only known 
gauge invariant operator to extract an interquark potential. 

At finite temperature, Euclidean time is compactified to the interval 1/T = aN t , 
with periodic boundary conditions for the link variables. The choice to = l,t = N t 
in Eq. (|5]) turns L(x) into a Polyakov loop encircling the torus in the time direction, 
transforming as £ 9 (x) = (7(x)L(x)^ _1 (x). In this case, the only gauge invariant object 
is the correlation of TrL, which is related to the color averaged potential, i.e. the 
weighted sum of the singlet and adjoint channels || ||, 

e -N tVav (r) = Tr ((P 1 + P ad )Q g )/Tr(P 1 + P od ) 

1 1 N 2 — 1 

= ^(TrLt(x)TrL(y)) = ^e"^« + £L^-^\ (9) 

However, temperature can be made arbitrarily small by increasing N t at fixed lattice 
spacing. The fact that there are gauge invariant operators for the singlet and the 



3 



average potential then suggests that there should also be a gauge invariant formulation 
for the adjoint channel. 

The problem can be solved by dressing the source fields with gluon "clouds" , whose 
fields transform in a fundamental representation || [Tj. Construction of such fields is 
not unique. One possibility is to use the eigenfunctions of the spatial Laplacian defined 
on every timeslice, 

-{^[U]) ap f^(x) = X n f^(x), A">0. (10) 

The latter is a hermitian operator with a strictly positive spectrum, whose eigenvectors 
have the transformation property f( n > 9 (x) = g(x)f( n \x), and may be combined into a 
matrix Q(x) G SU(N) following the algorithm given in ||10|| . It is a non-local functional 
in the sense that it depends on all links in a given timeslice. Since Eq. ( |10"D only deter- 
mines eigenvectors up to a phase, this leaves a remaining freedom in Q(x). In the case of 
£77(2), all eigenvalues are two-fold degenerate due to charge conjugation, and the two 
vectors to the lowest eigenvalue are combined into Q, which then is determined up to a 
global SU (2) rotation h. For £77(3) there is no degeneracy of the eigenvalues in general. 
In this case one solves for the three lowest eigenvectors to construct the matrix Q, which 
is then determined up to a factor h = diag(exp(ko>i), exp(icj 2 ), exp(icj 3 )), Yk u i = 0- 
This may be summarized by the transformation behaviour 

Q 9 {x) = g(x)Vt(x)h\t), (11) 

where h(t) is free and may be different in every timeslice. 

Let us then consider a dressed static quark by coupling it to such a gluonic field, 
ip n (x.) = f2^(x)'0(x). Replacing ip(x.) —* ip n (x) in Eq. ([|), one observes the desired 
transformation behaviour for all x, y, 

(P^(x)^(y)) 9 = (P^(x)^(y)) 

(P a ^ n (x)^(y)) 9 = -2£~(il> a t(x)htT*hlP(yj). (12) 

In the first line we have a singlet state. The fields Q(x) now play the same role as 
the spatial lines in the Wilson loop, allowing for flux between the sources to com- 
pensate their transformations. However, with this construction we have in addition a 
locally gauge invariant adjoint state, which transforms as an (N 2 — l)-plet under the 
residual global symmetry of the Laplacian. Time propagation of the dressed charge is 
correspondingly described by the dressed Wilson line 

L(x) = ft f (x, t )£(x)fi(x, t), (13) 
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which gauge transforms as L 9 {x) = h(t)L(x)hl(t). The correlator (TrLt(x)L(y)) is 
gauge invariant, and hence constitutes an alternative expression to extract the singlet 
potential. By means of the transfer matrix formalism JTT], [T2| the Euclidean expecta- 
tion values may be converted to traces over quantum mechanical states in a Hilbert 
space formulation, elucidating the quantum mechanical interpretation of the correlation 
functions. Comparing Wilson loop and the correlator of dressed lines, one has 

(W(r,t)) = Z- l fr{f N ^ a (x)U aP (x,y)My)F 
(TrZt(x)Z(y)) = 

Z~ l Tr {f JV -^ Q (x)^ 7 (x)^ /3 (y)^(y) f « $(y)0L(y)^(x)^(x)}. (14) 

where tp, $ are annihilation and creation operators for a static quark, (l acts as a 
multiplication operator and T denotes the transfer matrix. Inserting complete sets of 
eigenstates, it is evident that both expressions decay exponentially with the spectrum 
of the Hamiltonian, yielding identically the same singlet potential as their ground state. 

Next, consider the finite temperature case, where Eq. fljjD becomes a dressed Polyakov 
loop. Since TrL(x) = TrL(x), the expectation value (TrL) remains an order param- 
eter for the deconfinement transition of the pure gauge theory^ Moreover, the average 
potential is identically the same as that defined from L. Substituting dressed Polyakov 
loops for undressed ones in Eqs. (|7|), we thus obtain gauge invariant, non-perturbative 
definitions for the singlet, adjoint and average channels separately. 

How do the results depend on the particular construction of the functional f2[£7]? 
Since Q G SU(N), the dressing of the source may also be viewed as a gauge trans- 
formation by $l>(x), and the corresponding operator may be written as if) = Cvif> = 
R~ 1 (Q)ipR(Q), where R(g) is a unitary operator implementing a gauge transformation 
by g(x). The Hilbert space expression for the singlet potential then takes the form 

(TrLt(x)L(y)) = Z^T r {T^ a (x)^(y)i?(0) J* iT^^y^x)}. (15) 

The above expression is identical to the undressed one up to a similarity transforma- 
tion of the transfer matrix, which preserves the spectrum as well as the norm of the 
eigenstates. (Note that R(Q) does not commute with T because Q[U] depends on the 
link variables). Hence, the singlet potential can in principle be extracted by means of 
any Q[U], that is local in time and transforming as in Eq. ([11]). 

According to the above, the dressed Wilson line may equivalently be viewed as a 
Wilson line brought to Laplacian Coulomb gauge by a gauge transformation g{x) = 

2 Note that even the untraced loop £(x) retains the same transformation behaviour under the 
center of the gauge group as £(x), since the functionals fl(x) are solutions of the spatial Laplacian. 
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fl^x). In the language of gauge fixing the previous statements may then be rephrased 
as follows: the gauge fixed Polyakov loop correlation function in the singlet channel 
falls off with gauge invariant eigenvalues of the Hamiltonian. These may be extracted 
in any unique gauge that is local in time, i.e. preserves the spectrum of the transfer 
matrix. This is true, for example, for the Coulomb gauge, but not for the Landau or 
complete Laplacian gauge, which depend on all time-like links as well. 

In the remainder numerical results will be presented supporting these theoretical 
statements and illustrating their numerical feasibility. This first exploratory study 
is done for SU(2) in 2+1 dimensions, for its significantly lower numerical cost and 
fast continuum approach. In this case the coupling constant g 2 has dimension of 
mass and provides a scale. As a consequence, the short distance Coulomb part of the 
potential defined through the Wilson loop is logarithmic, and in next-to-leading order 
a contribution to linear confinement is obtained even in perturbation theory, 



O r ) = i^V) + a r + (g 6 r% (16) 

Z7T 



where (J pert however falls short of the non-perturbative string tension by a factor of 
about 2/3 fpL3| . Since the group theory is unaffected by the number of space dimensions, 
the leading order perturbative relation between singlet and adjoint Coulomb terms, 
Eq. (fj), carries over to this case (where the SU(2) adjoint channel corresponds to a 
triplet, V a d = V 3 ). The lattice gauge coupling is given by (3 = 4/(ag 2 ). 

First, it is instructive to test the numerical feasibility of the dressed Wilson line 
by comparing its zero temperature correlation function with the Wilson loop, as in 
Eq. (|H]). Fig. [I] shows the singlet potential extracted from those correlation functions, 
and numerically confirms their analytically established equality. (Numerically this has 
also been observed for Wilson loop and lines in the adjoint representation 0]). 

Next, we consider finite temperatures and dressed Polyakov loop correlations. The 
deconfinement transition of the 2+1 dimensional gauge theory has been studied in 
15|,||16||. The relationship between the inverse critical temperature in lattice units, 



N£ = l/(aT c ), and the lattice gauge coupling j3 is well fitted by N£(/3) = (/? — 0.37)/1.55 
171 . For the following calculation of the static potentials we work at (3 = 9, for which 



the deconfinement transition is between iV t c = 5, 6. The three potentials are shown for 
decreasing temperatures in Fig. 

The results clearly confirm the perturbative expectation that the triplet potential is 
repulsive at short distances. In the deconfined phase all potentials assume a Yukawa 
form and merge at a common constant value at large distance, i.e. all color interac- 
tions are screened. In the confinement phase, however, the triplet potential becomes 
attractive at larger distances, while staying always above the singlet potential for the 
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Figure 1: Singlet potential extracted from the Wilson loop (W(r,t)) (open squares), and 
correlations of dressed Wilson lines (L(x)L(y)) (crosses), for (3 = 9, L 2 ■ N t = 32 3 . 



distances covered in the simulation. This is in accord with the fact that no stable 
charged states are observed empirically. Nevertheless, at low temperatures a potential 
well seems to form in which a metastable bound state appears possible. An investiga- 
tion of this phenomenon in 4d SU(3) would be particularly interesting for heavy quark 
phenomenology, such as J/\l/-production through an octet channel. A nice consistency 
check for the calculation is provided by comparing the singlet potential in the low tem- 
perature case with the zero temperature potential obtained from the Wilson loop. The 
latter is indicated in the figure for T = 0.35T C by the dotted line, which is approached 
by the low T Polyakov loop data from above, as expected on physical grounds. 

As reported in the non-locality of the functionals Q[U] has them project predom- 
inantly on flux tubes encircling the periodic spatial volume, hampering the observation 
of localized states by large finite volume effects. Such effects are absent in the cur- 
rent calculations, as is demonstrated in Fig. |^ in the low T plot, where the data for 
L = 32, 48 yield entirely consistent results. This is easily understood recalling that the 
quark- ant iquark system is spatially extended, and the non-locality of the Q(x) actu- 
ally enables flux tubes to connect the pair over large distances. The flux represented 
by Q(x) also closes through the periodic boundary conditions. As a consequence, we 
really measure a superposition of exp —N t V(r) and exp — N t V(L — r). However, for 
our choices of N t , L, the second term is exponentially suppressed compared to the first 
one, and no finite volume effects are observed. 

In order to illustrate the independence of the calculation of the particular choice 
of Q, the T = 0.7T C plot in Fig. ^| shows an additional calculation, where Q[U] is 
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Figure 2: Singlet, triplet and average potentials calculated for different temperatures at 
(3 = 9. At T = 0.7T C , the crosses denote data obtained in standard Coulomb gauge. At 
T = 0.35T C; filled symbols denote data obtained on a 48 2 • 16 lattice. The dotted line 
marks the T = Wilson loop data from Fig. [J. 



the functional to fix standard Coulomb gauge, i.e. it minimizes the function R[U] = 
X)a;i[l — 1/iVTr Q> (x)Ui(x)Q(x + i)] in every timeslice. The corresponding data points, 
shown as crosses in the plot, coincide with those from the previous calculation. The 
reader is reminded that the lattice Coulomb gauge in principle is afflicted by Gribov 
copies. Since we know from our analytical considerations that the potentials have to 
be independent of the choice of Q, this numerical result may be interpreted as a test 
for the effect of Gribov copies in the standard Coulomb gauge. 

Having separate non-perturbative results for the singlet and adjoint channel poten- 
tials, they should approach the leading order perturbative prediction, Eq. (g), in a 
regime where perturbation theory is valid. This is the case at short distances, and ex- 
pected to extend to larger distances in the deconfmed phase. Since the static potential 
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includes an unphysical divergence which is different in continuum and lattice regular- 
izations, we compare the physical forces F t = dVi/dr. Fig. [3] shows results obtained in 
the deconfined phase on a finer lattice, by employing the discretized forward derivative 
Fi(r) = A r Vi(r) = (Vi(r + a) — V(r))/a. Satisfactory agreement with the leading order 
perturbative result is observed. 



x aF, 
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Figure 3: Comparison of the forces in the singlet and adjoint channels, computed at 
(3 = 18, L 2 ■ N t = 48 2 • 8, corresponding to T = 1.42T C . 

To conclude, in following the methods of [§, [7], we have shown that a dressed static 
quark- ant iquark pair can exist in locally gauge invariant singlet and adjoint states, the 
latter transforming as a global (N 2 — l)-plet. Calculation of the correspondingly dressed 
Polyakov loop correlation functions then permits the non-perturbative extraction of the 
singlet, adjoint and average potentials which are each eigenvalues of the Hamiltonian 
and hence gauge invariant. Numerical results for 2+1 d SU(2) gauge theory provide a 
consistent picture, in agreement with perturbation theory at small distances and large 
temperatures. Moreover, the calculation allows to observe how the finite temperature 
singlet potential extracted from the dressed Polyakov loop approaches the Wilson loop 
results in the zero temperature limit. The adjoint channel potential is repulsive at short 
distances. At zero temperature it becomes confining at large distances, allowing for a 
metastable bound state which should be of interest for heavy quark phenomenology. 
The considerations in this work can be applied to construct gauge invariant symmet- 
ric and antisymmetric quark-quark potentials as well. Apart from their relevance for 
heavy quark physics and the confinement problem, our results should also offer new 
ways to study the forces and excitations in hot and dense QCD. 
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